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Progress in automated perturbation theory for heavy quark physics 

Matthew A. Nobes a , Howard D. Trottier a 

a Physics Department, Simon Fraser University, Burnaby, BC, Canada 

We review our perturbative techniques for improved heavy quark actions. A new procedure for computing 
improvement coefficients is suggested, where the continuum limit of a lattice-regularized theory provides the 
matching conditions. We also use a gauge-invariant infrared regulator that is well suited to higher-order calculations 
in lattice gauge theory. We report on preliminary tests of our method, as well as a possible way to reduce systematic 
errors in these calculations. 
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1. Introduction 

Improved actions are widely used in heavy 
quark physics. Examples include NRQCD pQ 
which is very successful in describing bb systems, 
and the Fermilab action [2] which has been ap- 
plied to both b and c quark systems. With the 
MILC results for unquenched light quarks, preci- 
sion determinations of important physical quanti- 
ties (i.e. fs) are possible with these heavy quark 
actions 

One of the dominant systematic errors that re- 
mains with improved actions is the determina- 
tion of the improvement coefficients. The hyper- 
fine splittings of the J/ip system illustrates this 
clearly. While there have been good results from 
heavy quark simulations, these hyperhne split- 
tings remain significantly underestimated with 
tree-level improved actions. This underestimate 
is seen with both the Fermilab and NRQCD ap- 
proaches to heavy quark physics. Unquenching 
(010) does not seem to resolve it (the Fermilab 
determination is still w 25% low). We expect that 
this underestimate should be greatly reduced by 
perturbative matching of the S • B operator and 
the use of a highly improved action j^j. 

In order to determine this action coefficient 
(cb) we adopt the simplest matching strategy 
possible. We will compute the amplitude for 
quark scattering off of a background chromo- 
magnetic field in lattice perturbation theory, and, 
to the same order in perturbation theory compute 
the continuum result, using the same IR regula- 



tor. The matching is simply a matter of subtract- 
ing one from the other, and tuning cb to insure 
that the difference vanishes to whatever order is 
desired 

In principle, this procedure should work to 
whatever order in perturbation theory we wish. 
In this report we will discuss the matching to one- 
loop order, however for true high precision results 
the matching should be done to two loops. This 
is because at typical lattice spacings ay (1/°) ~ 
0.2. Our ultimate goal is high precision results, 
so we want to organize our calculation in a way 
which makes going to two loop order as easy as 
possible. 

We have discussed our approach to lattice per- 
turbation theory at length in prior LATTICE pre- 
sentations ( , |E] > EH ) • The core of our technique 
is to generate Feynman rules automatically using 
the Luscher - Weisz algorithm ^U]. The primary 
advantage of this method is that we can change 
the form of the action easily, despite the complex- 
ity of the rules. 

Many of the individual diagrams we are calcu- 
lating are infrared divergent, to regulate this we 
use twisted boundary conditions. This has the ef- 
fect of turning integrals over the three momentum 
into restricted sums, as in the following example: 
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Here Xk is a veto function which excludes many 
lattice momenta, including the IR zero modes. 
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2. Lattice-to-Lattice Matching 

By design, the lattice and continuum theories 
should be identical in the infrared. The lattice 
matching computations that have been done to 
date have typically used a gluon mass to regulate 
the IR in both the lattice and continuum theories. 
This method could work for our calculation, how- 
ever, it may pose significant problems with gauge 
invariance at higher order. 

One solution is to use twisted boundary con- 
ditions, which provide a gauge invariant IR reg- 
ulator. However this regulator is somewhat dif- 
ficult to combine with dimensional regularization 
for continuum calculations. This can be avoided 
by rethinking the continuum UV regulator that 
we use. We are calculating matching coefficients, 
so there is no particular reason to use dimensional 
regularization. In fact, we can use any UV regu- 
lator we want. 

Following a suggestion of Peter Lepage we have 
simply used a lattice cutoff to regulate the con- 
tinuum theory. This allows us to do the "contin- 
uum" side of the computation using the same lat- 
tice techniques as we have developed for the heavy 



Figure 2. Wavefunction renormalization for light 
Wilson and SW quarks 



quark actions. This procedure is straightforward, 
we can pick a simple lattice theory (Wilson glue 
+ naive quarks) and run the lattice spacing down. 
As long as the spacing is made small enough, this 
theory will be very close to the continuum theory. 
Then we can subtract off our heavy quark results, 
and get the matching coefficients. This lattice-to- 
latticc technique should work for matching calcu- 
lations within QCD itself, such as the perturba- 
tive coefficients for operators in the lattice action. 
This method does not apply for example to ma- 
trix elements of operators in the effective Hamil- 
tonian for the weak interactions. 

As we discuss below, for one loop matching 
computations this naive procedure works. To two 
(and higher) loops the "continuum" and lattice 
theories will both have to be renormalized to give 
the same physical values of the various param- 
eters, however this should be a straightforward 
application of renormalized perturbation theory. 
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3. Testing 

To test this procedure we do the computation 
of c^ b m the static quark limit. This calcula- 
tion has previously been performed, using a gluon 
mass as the IR regulator, and matching to the 
dimensionally regularized continuum theory, by 
Flynn and Hill JT]. We compute the amplitude 
for quark scattering off of a background field for 
both the heavy quark and "continuum" theories. 

In the "continuum" case, we're using units in 
which a = 1 so we want to run all other scales in 
the problem down to zero. For matching, we use 
naive quarks, so the differences from the contin- 
uum limit go like 0(m 2 ). 

Figure 2] shows the results of this lattice to 
lattice matching for various naive quark masses. 
The black squares are the results of ^Ij • Clearly, 
in the continuum limit, L — > oo and m — » 0, we 
reproduce the Flynn-Hill result. Figure ^ shows 
that our matching procedure works, and that the 
errors will be small enough to reliably extract the 
dependence on the ultraviolet cutoff. 

In addition to this test, we have throughly 
tested our implementation of the Fcrmilab 
fermion action. To do this we have duplicated all 
of the major results in |12| for the one loop rest 
mass and wavefunction renormalization. Figure 
Elshows the (one loop) wavefunction renormaliza- 
tions of Wilson and clover quarks, at fixed lattice 
size. The logarithmic dependence on Mq agrees 
with |T2]. 

4. Periodic Boundary Conditions 

We have investigated other means for control- 
ling infrared divergences. One possibility is to 
integrate the difference between the lattice and 
the "continuum" expressions for a matching co- 
efficient. This difference is infrared finite, and 
hence could be computed without any infrared 
regulator. 

To test this method, we have calculated the dif- 
ference in fermion wave function renormalizations 
^Wilson _ %sw Qver a ran g e f fermion masses, 

without any infrared regulator. In the mo — > 
limit we obtain 0.228(1) which is in reasonable 
agreement with 0.023083(7) which is the result 



presented in ^2]- We are currently investigating 
this technique for determining cb- 

5. Conclusions 

Doing high precision QCD with improved 
heavy quark actions requires two loop determi- 
nations of the improvement coefficients. In par- 
ticular the coefficient of the a ■ B is needed to 
resolve the underestimate of the hyperfine split- 
ting we have discussed. 

To get to two loops it is important to orga- 
nize the calculation in as efficient a manner as 
possible. To this end we have introduced a pro- 
cedure to match two lattice theories onto each 
other, where the continuum limit of one lattice 
regularized theory provides the matching crite- 
ria for the improvement of the other. Addition- 
ally we have used twisted boundary conditions to 
provide a gauge invariant regulator. As a prelim- 
inary test of these techniques we have reproduced 
an earlier result for the matching of a ■ B in the 
lattice HQET theory. We are currently applying 
these methods to improvement of the Fermilab 
and NRQCD actions. 
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